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This work deals with bifurcation and pattern changing in models described by two real scalar 
fields. We consider generic models with quartic potentials and show that the number of independent 
polynomial coefficients affecting the ratios between the various domain wall tensions can be reduced 
to 4 if the model has a superpotential. We then study specific one-parameter families of models and 
compute the wall tensions associated with both BPS and non-BPS sectors. We show how bifurcation 
can be associated to modification of the patterns of domain wall networks and illustrate this with 
some examples which may be relevant to describe realistic situations of current interest in high 
energy physics. In particular, we discuss a simple solution to the cosmological domain wall problem. 



I. INTRODUCTION 

Bifurcation is a well known phenomenon, which occurs 
in several distinct areas of nonlinear science [1 . Usually, 
it can be controlled by some forcing parameter, which 
drives the system to change from a given behavior to an- 
other. Interesting examples can be found, for instance, 
in a competition between the roll and hexagon structures 
which appear in an experiment of thermal convection 
[2], in nematic liquid crystals subject to rotating mag- 
netic fields where spiral waves were observed [3], and in 
periodically forced nonlinear oscillatory systems, where 
patterns are seen to spring induced by bifurcation, as 
reported for instance in [H |5] . 

A distinguishing feature of these studies is that they 
are all non-relativistic. However, bifurcation also appears 
in various high energy physics models where static do- 
main wall solutions are shown to bifurcate under the vari- 
ation of model parameters [6 . Other examples are given 
in [7], relevant for the cosmological evolution in brane- 
world scenarios, in [8j, associated with critical gravita- 
tional collapse, and in [9l [T^, in the context of hybrid 
inflation. 

In the present work, we take a specific direction by 
studying bifurcation in a class of relativistic models de- 
scribed by two real scalar fields, (j) and x- The basic moti- 
vation is to study the possibility that a network of domain 
walls may bifurcate into different configurations, giving 
rise to modifications of the network pattern. If we choose 
two minima arbitrarily, it may be possible that they are 
connected with two or more distinct orbits. When this 
happens the system develops a bifurcation, since one can 
go from a given minimum to another one, following two 
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or more distinct paths. In one dimension, given three 
minima j, /c, and /, the direct jl path is energetically 
favored over the path jk + kl if the energy density or 
tension of the {jl) structure is lower than the sum of the 
other two tensions, that is rji < rjk -\- Tki. A bifurcation 
occurs when the inequality becomes an equality, since we 
can go from j to / following the direct jl path, or alter- 
natively visiting k through the path jk -\- kl^ at no extra 
energy cost. In more than one dimension the analysis 
needs to take into consideration the intersection angle of 
the domain walls. 

In this paper we consider generic models with quartic 
potentials with or without a superpotential. In the for- 
mer case the analysis is simpler since the tension of the 
domain walls associated with a given BPS sector is fully 
determined by the values of the superpotential, x)^ 
at the minima of the potential. We take x) be 

a polynomial function of degree 3 which gives rise to a 
quartic potential, x) (note that for superpotentials 
of order lower than 3 there are no domain wall solutions) . 

The situation is more complicated when the two-field 
model is not the bosonic portion of some supersymmet- 
ric theory. In this case, there is no superpotential and 
the determination of the domain wall tensions is not so 
straightforward. However, an interesting procedure was 
followed in pTj where a model (the so-called BBL model) 
with two real scalar fields, (j) and x, was constructed un- 
der the assumption that its potential could be written as 
the sum of two quartic potentials, V{(j)) and V{x)^ plus 
a small interaction term. Although the model as a whole 
has no superpotential, the importance of the interaction 
term can be controlled by adjusting a small parameter, 
e. The domain wall tensions could then be determined 
analytically for \e\ <C 1. 

The evolution of the macroscopic properties of domain 
wall networks is strongly dependent on whether or not 
junctions are present and, in the former case, on the di- 
mensionality of the junctions. Although the presence of 
junctions is not enough to make domain walls a viable 
dark energy candidate, their impact on the dynamics of 
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the network needs to be taken into account in any quan- 
titative study of cosmological consequences of domain 
wall evolution [12]. In this paper we perform domain 
wall network simulations of a new class of one-parameter 
multi-tension domain wall models, showing that it has a 
number of interesting properties. In particular, we dis- 
cuss a simple solution to the cosmological domain wall 
problem which arises naturally in this class of models. 

In this paper we study bifurcation and pattern chang- 
ing in models described by two real scalar fields, with or 
without a superpotential. The outline of this paper is as 
follows. We start by reviewing the BBL model in Section 
II, analytically and numerically calculating the domain 
wall tensions and comparing our results with analytical 
results, valid for \e\ <C 1, obtained in [TP (see also [13 ). 
We discuss how bifurcation arises in this model as well 
as the corresponding implications for the dimensionality 
of the junctions which may be present when we consider 
more than one spatial dimension. In section III we study 
the other possibility, that is, we consider models which 
support bifurcation and are described by a superpoten- 
tial. We investigate this case in detail because it is ex- 
pected to be of more general interest, since it can be seen 
as the bosonic portion of a larger, supersymmetric the- 
ory. We also discuss some implications of our findings in 
the light of the cosmological domain wall problem. We 
end the paper in section IV with some comments and 
conclusions. 



II. ABSENCE OF SUPERPOTENTIAL 

Let us start by reviewing the two field model intro- 
duced in [11 . Consider two real scalar fields (j) and x 
with the Lagrangian 
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-^:^d,xd^X-V{^,x). (1) 



where the potential V is given 

V{cP,x)=V{cP) + V{x) + eu{^,x), (2) 

with e being a small parameter, which control the inter- 
action between the two fields. 
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and r is a real number. The numerical factors are intro- 
duced for convenience. We are considering dimensionless 
variables, that is, we have rescaled both the fields and 
the space and time coordinates appropriately such that 
they are now dimensionless quantities. 

If £ = then the fields, (j) and x^ coupled and 

the minima of the potential are at = ((^i, Xi) with = 
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FIG. 1: Comparison of the values of Te determined analyt- 
ically, using either the straight line approximation (curved 
solid line) or the linear approximation (straight dot-dashed 
line), with the numerical solution (filled circles), as a func- 
tion of e. The curved dashed line represents the value oiTd/2 
calculated using Eqn. (6b). We have taken r — 1. 



{—r^r)^ V2 = {—r^—r)^ = (r, — r) and = (r, r). For 
e very small the minima move slightly to new positions 
vf. In this limit the tensions can be calculated as [M] 



where 
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(5b) 
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8r^/3 are the tensions for the 



edge and the diagonal solutions (j)^ and (/)^, respectively 
(with similar expressions for the x-field) . The exact solu- 
tions with 8 = can thus be used to obtain the tensions 
r| = 4r2(l + 7£rV4)/3 and r| = 8r2(l + er'^)/3 (see 
Refs. J.3 for more details). 

By considering straight line orbits connecting two dis- 
tinct minima at the edge and diagonal sectors, we get the 
corresponding tensions 



(6a) 
(6b) 
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In fact, due to the symmetry of the scalar field potential, 
the straight diagonal orbit is a solution of the equations of 
(4) motion, so that Td = rf corresponds to an exact solution. 
Although the edge solutions are not straight lines, Eqn. 
(6a) provides an approximation to the real tension up to 



first order in e. We have also computed Tg by solving the 
equations of motion for (j) and x numerically (see Ref. 
p!2l for details about the numerical code). 

In Fig. 1 we compare the values of Te determined nu- 
merically as a function of e (filled circles), with the an- 
alytical expectations for \e\ <C 1 determined using either 
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FIG. 2: The left panel illustrates that for e < the model 
forms stable X-type junctions. The right panel shows that the 
formation of y-type junctions is possible for s > 0. We took 
€ — —0.2 and e — 0.2 in the numerical simulations illustrated 
in the left and right panels respectively. 



the straight line approximation (Eqn. ( |6a[ ) - curved solid 
line) or the linear approximation (Eqn. ( [5a| ) - straight 
dot-dashed line). The curved dashed line represents the 
value of Td/2 calculated using Eqn. (|6b|. Here we made 
the choice r = 1. We see that the analytical results for 
Tg are a good approximation of the exact solution for 
small values of becoming significantly less accurate 
for \£\ > 0.2. We also see that Td/T^ is greater or smaller 
than 2 depending on whether e is smaller or greater than 
zero, respectively. To investigate bifurcation, in one spa- 
tial dimension, we have to compare the diagonal tension 
with the sum of two edge tensions, and so we see that 
for £ < the model relaxes to the edge solutions, and for 
s > it relaxes to the diagonal configuration. Thus, the 
sign of £ induces bifurcation. 

The above discussion in one dimension is crucial to de- 
termine whether F-type or X-type junctions are favored 
in more than one dimension. If r^^/rg > 2 {rd/re < a/2) 
then only X-type (F-type junctions) will be formed. 
There is an intermediate case corresponding to a/2 < 
'?"d/'?"e < 2 in which both stable F-type and X-type junc- 
tions are possible (depending on the intersection angle 
between the walls). However, if one starts with a X- 
type junctions where four walls intersect at right an- 
gles, it will be stable for Td/re > a/2 and unstable for 
^d/'^e < a/2- Let us consider the intersection angles Oi 
between the walls that surround each domain with min- 
imum Vi {i = 1,2,3 or 4). Notice that 6i = 6s and 
O2 = O4. For e < (X-type junctions) equilibrium is 
possible even if Oi is different from O2. Otherwise the do- 
main walls would intersect at right angles which does not 
need to be the case, as it is clearly shown in the left panel 
of Fig. 2. For e > stable X-type junctions can still ex- 
ist if 2reCos(^2/2) < Td (assuming that O2 < Oi). In 
this case, O2 = 7r/2 is the maximum possible value of the 
intersection angle. Consequently, we see that stable X- 
type junctions are also possible for any a/2 < Td/re < 2, 
as long as O2 > 2 ai cos {rd/{2re)) . For smaller values of 
^2, 1^-type junctions will form. 

We illustrate some of these points in Fig. 2 where we 
plot the results of two different simulations of the BBL 
model in two spatial dimensions. We take e = —0.2 and 



e = 0.2 in the left and right panels, respectively. In the 
first case, Td/re > 2 and consequently F-type junctions 
are unstable which leads to the formation of X-type junc- 
tions. In the second case, Td/Te < 2 and consequently the 
formation of F-type is no longer forbidden. The mul- 
tiplicity of the junctions is one of the most important 
variables affecting domain wall network evolution and it 
is an essential aspect to consider in any quantitative as- 
sessment of the cosmological implications of domain wall 
networks [T2j. 



III. PRESENCE OF SUPERPOTENTIAL 

In the present work we deal with bifurcation and pat- 
tern formation in two-field models, characterized by the 
real scalar fields (j) and described by potentials con- 
taining up to fourth-order power terms. 



A. Generic models with cubic superpotentials 

We shall now consider a model with two real fields 
characterized by the potential 



where W = W{<p,x) is the superpotential and 

dW dW 

d(f> dx 



(7) 



(8) 



The equations of motion for static solutions have the gen- 
eral structure 



<fx 



= W4,W4,4, + W^W^4,, (9a) 
W^W^^ + W^W^^ . (9b) 
They are solved by solutions of the first-order equations 



dx ^ 



(10a) 
(10b) 



In the following we shall drop the ± sign. It will be suf- 
ficient to realize that for each solution of the first order 
equations, (^(x),x(^)), there will also be another solu- 
tion given by ((/)(— x), x(—^))- The minima of the po- 
tential are the points in field space (^, x) which solve the 
two algebraic equations VK^ = and = 0, and so they 
are critical points of W. We label the minima as j, /c, / , 
with j = {/^j^Xj)-) 6tc. Each pair of distinct minima j and 
k {k 7^ j) forms a topological sector. A static solution 
is represented by (0(x),x(^))- If it solves the first-order 
equations and connects two distinct minima in field space 
then this solution is a minimum energy configuration in 
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this sector. It corresponds to topologically stable solu- 
tions, usually named BPS states in the supersymmetric 
model, corresponding to a topological BPS sector - see 
[TF for more details. A non-BPS sector is a topologi- 
cal sector which has no solution coming from solutions 
of the first-order equations. This is easy to identify be- 
cause the tension associated with a BPS sector is given 
by Tjk = \W{j) - W{k)\ 7^ 0, and so a non-BPS sec- 
tor is specified by the minima / and m, say, which obey 
W{1) = W{m). 

We may write a generic potential of order 2{d — 1) for 
the scalar fields (j) and x 



0<a+6<2((i-l) 



(11) 



The number of coefficients, Cab^ of the polynomial po- 
tential is given by 



Nd = 2d^ -d, 



(12) 



so that A/'2 = 6 and = 15. On the other hand, if 
the above potential can be obtained from a polynomial 
superpotential of order d. 



(13) 



l<a+h<d 



using Eqn. ([7|), then the number of independent coeffi- 
cients of the polynomial potential reduces to 



(14) 



that X2 = X4 = and take 1/2 = ct and = b. For 
a —b the most generic cubic superpotential satisfying 
all of the above conditions is 



4ab 



2x 



4ab 2 4 1 



a -\- b a -\- b 3 a -\- b 

If a = —b the superpotential is given by 

Wi^,x) = ^(-36^X + 36VX + X') 



(15) 
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+ 3^x^ + 6 



(16) 



The standard BNRT model fi6\ is obtained by making 
the choice Ai = and A2 = — 1 but other choices of Ai 
and A2 also do not modify the positions of the minima 
'^2 = (0, —b) and V4 = (0, b). 

In this paper we shall take Ai = 6^ and A2 = so that 
the superpotential becomes 



W{<P,x) = -bx- 



(17) 



with 6 > 0. This particular choice is equivalent to a ro- 
tation of the original BNRT model by 7r/2 followed by 
a rescaling of the fields by a constant factor and con- 
sequently does not have any further implications. This 
model has a Z4 symmetry if 6 = 1 which is broken into 
Z2 X Z2 for b^l. 

In this case the first order equations are given by 



Here, Cah are the coefficients of the polynomial super- 
potential. We see that N2 = ^ and A/'a = 9. How- 
ever, some of these are not real degrees of freedom in 
the sense that there are changes in the polynomial coefi- 
cients which do not affect the ratios between the domain 
wall tensions (assuming that the model has domain wall 
solutions). These ratios are invariant under a change of 
the overall normalization of the potential V ol^V a 
rescaling of the fields (0, x) ^ (<^2^, <^2X)7 a translation 
(^7 x) ~^ (0+0*7 X+X*) oi" a rotation by an angle d where 
a\ > 0, 0^2, 0* and x* are real numbers. It is clear from 
the discussion that the case with = 2 is not interesting 
because the lowest polynomial order of W for which one 
can obtain domain walls interpolating between different 
minima is d = 3. In this case the initial number of coeffi- 
cients, 9, is reduced to 4. In practice one can arbitrarily 
choose the positions of two of the minima at '^i = (— 1, 0) 
and i'3 = (1,0) as well as the normalization of the poten- 
tial and the remaining 4 coefficients are associated with 
the positions of the other two minima V2 = (^2,^2) and 

= (x4,?/4). 



B. A class of models with bifurcation 

For simplicity we shall now consider a class of models 
with the two minima, V2 and v^^ lying in the yy axis so 



dx 
dx 



^bcl^x. 



The superpotential at the four minima is given by 
W{vi) = 0, W{v2)='^b\ 
W{v3) = 0, W{v4) = -lb^. 



(18a) 
(18b) 

(19a) 
(19b) 



Thus, the model has one topological sector that is non- 
BPS, the sector vi ^ vs^ and five BPS sectors, one with 
tension r24 = 46^/3, and four with the same tension ri2 = 
Ti4 = r23 = r34 = 26^/3. Also, the BPS sector V2 ^ v/^ 
has an infinity of BPS solutions, degenerated to the same 
tension r24 = 46^/3. A vertical straight line BPS solution 
is given by 



(j){x) 



0, 



-6tanh(x) . 



(20a) 
(20b) 



The horizontal solution connecting the minima vi and ^3 
is non-BPS and is given by 



(j){x) = tanh(6x) , 

x(^) = 0, 



(21a) 
(21b) 
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FIG. 3: In the top panel we plot the BPS tensions as a func- 
tion of b. Solid, dashed and dot-dashed lines represent the 
tensions T24, n2, and ri4, respectively. In the lower panel we 
plot the sum of the BPS tensions ri4 + T34 and ri2 + r23 for 
6 > as well as, rfl, for the straight line non-BPS solution. 
Note that if 5 > then r24 = T12 + n4 and that T14 = T34 
and T12 = T23 for any value of b. Also, notice that rf^ T13 
if 6- 1. 



perpotential which is given by 

W{4>, x) = x(-b+l{l- b)x + b<p^ + ^x') , (22) 
SO that the potential is 

V{ct>, x) = \ {{2b4>xf + i-b + (1 - b)x + b<p^ + x^f) ■ 

(23) 

We see that if 6 = then is at the origin, that is V4 = 
(0,0), and the minima vi^ ^3, and V4 are aligned. When 
b = 1 the four minima are at the vertices of a square. 
Other special points are (see below): b = —1/3 giving 
V4 — (0,-1/3), b = —1 making V4 = (0,-1) coincident 
with and b = —3 giving V4 = (0, —3). 

In the case of a generic 6, the first-order equations are 

g = ^b<Px, (24a) 

g = -b+{l-b)x + b<l>^+x^. (24b) 

There is a particular solution with = 0, which connects 
the minima V2 and ^4, for which we have 

|^(.+ i(l-'.))'-i(l + '.)^ (25) 

which has the explicit solution 

X{x) = -1(1 - 6) - 1(1 + b) tanh((l + b)x/2) . (26) 

Note that in the b ^ —1 limit the minimum V4 collapses 
into eliminating the topological sector V2 ^ V4. 
The superpotential evaluated at the four minima gives 



with tension rfl = 46/3. This model has many interest- 
ing features. In particular, all the orbits can be found 
analytically in field space, as it was shown in Ref. [17]. 
Another feature of this model is that it has r24 = Ti2+ri4, 
independently of the value of b. This has been studied 
in detail in Refs. [151 IISl Hi] we shall now turn our 
attention to another one-parameter class of multi-tension 
domain wall models which has not appeared before in the 
literature. Unlike in the BNRT model, the non-diagonal 
sectors do not all have the same tension. We will see later 
that this property will be crucial in the light of a possible 
solution to the cosmological domain wall problem. 

The class of models we introduce in this paper is char- 
acterized by a cubic superpotential with minima at the 
values vi = (—1,0), V2 = (0,-1), '^3 = (1,0) and 
'U4 = (0,6). For b = 1 the model has a Z4 symmetry 
which is broken into Z2 for 6 7^ 1. 

Apart from an arbitrary normalization, the positions 
of the minima completely determine the form of the su- 



W{vi) = 0, W{v2) = -{l^3b) , 





(27a) 



W{vs) = 0, W{V4) 



6^(3 + 6), (27b) 



and for 6 = —1 we get that W{v2) = W{v4), in accor- 
dance with the above discussion. However, for any value 
of 6 7^ —1, the topological sector V2 ^ V4 is always a 
BPS sector, and the above particular solution (26) is al- 



ways a BPS solution. With (27), we can then collect the 



tensions in all the BPS sectors. They are given by 



ri2 

n4 



T23 

> 

T34 



1 
6 

1|3 

1 
6 



1^3b\. 



b'\3^b\. 



(28a) 
(28b) 
(28c) 



We plot the BPS tensions in the top panel of Fig. 3, to 
see how they behave as a function of b. Solid, dashed and 
dot-dashed lines represent the tensions r24, ri2, and ri4. 
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respectively. In particular, we note that r24 = Ti2+ri4 for 
b < —3 and b > —1/3, ri4 = ri2+r24 for —l<b< —1/3, 

'Tl2 = '?"24 + ^14 for —3 < 6 < — 1. 

For 6 = the three minima vi, ^3, and 'U4 get aligned, 
with W{vi) — ^^('^3) = W{v4) =0, so they cannot form 
BPS sectors. This case is uninteresting because we can 
make the (j) field constant, with the x fi^ld obeying 



such that for 6 = the superpotential has the form 



(29) 



(30) 



The b — 1 case leads to a similar situation, and the 
fields (j)-^ = (j)-\-x and (j)- = (j) — X decouple, but now they 
both have nontrivial behavior. When b = —1/3, we get 
that W{vi) = W{v2) = W{v3) = and ^(^4) = -4/81. 
In this case there are three BPS sectors with the same 
tension, ri4 = r24 = T34 = 4/81. For b = —3 we get 
W{vi) = W{v3) = W{v4) = and W{v2) = 4/3. Hence, 
in this case there are also three BPS sectors with the 
same tension, ri2 = r23 = r24 = 4/3. This is similar to 
the former case, for b = —1/3, with the minima V2 and 
V4 interchanged. 

In the sector vi ^ V2 (and similarly in the sector V2 ^ 
vs) we can find a straight line orbit connecting these two 
minima. In the sector V2 ^ V3 if we require that ^ = x+1 
then the potential becomes (see Eqn. (|23|)) 



1 



2_2 



v{x) = ^p\x^irx 

with p = (5572 + 6 + 1/2)1/2. rj.^^ solution is 



tanh {P x) , 



(31) 



(32) 



SO that the corresponding tension of the straight line so- 
lution is 



f3 



(33) 



There are similar expressions in the sector vi ^ V2 
(ti2 = T23). Interestingly, we note that, except for b = 
— 1/3, these two sectors are BPS sectors and they have 



the tension r23 < rl^. We compare this with Eqn. (|28a) 



to conclude that, in this sector, there are no straight line 
solutions to Eqns. (24a) and (24b), except when 6 = 1. 



Still, Eqn. ( 33 ) provides an upper bound to the tension of 
a static domain wall in the sectors vi ^ V2 and V2 ^ v^. 
In the case of BPS solutions this turns out not to be very 
helpful since the exact values of the domain wall tensions 
can easily be calculated. However, this is no longer true 
for non-BPS solutions and consequently we may use the 
above procedure in order to obtain useful upper bounds 
on the tension of the domain walls. 

For 6 7^ the sector vi ^ vs is always non-BPS. Let 
us find a straight line orbit connecting the two minima 




FIG. 4: Plots of the minima of the potential. The drawings 
start at the left top, with b given respectively by 6 = 1, 6 = 0, 
b — —1/3, b — —1 and b — —3. The colors and shapes of 
the points are used to identify distinct values of the superpo- 
tential at the corresponding minima. Solid and dashed lines 
indicate BPS and non-BPS sectors, respectively. For 6 = 0, 
the horizontal red and green lines show the continuum of min- 
ima with distinct values of the superpotential. 



vi and ^3 for 6 7^ 0. We consider = Q to get to the 
potential (see Eqn. (23)) 



(34) 



This means that the straight line non-BPS solution is 



(j){x) = tanh(6x) . 



with tension 



'13 



:4|6|/3. 



(35) 



(36) 



In the lower panel of Fig. 3 we plot the values of the 
sum of the BPS tensions ri4+r34 and ri2+r23 for 6 > as 
well as rfg, the tension corresponding to the straight line 
non-BPS solution. Note that the straight line solution 
is expected to become increasingly accurate in the limit 
6 ^ 1 (so that rfg ~ ri3 for b ^ 1). 

To visualize the behavior of this family of models, in 
Figs. 4 and 5 we depict the minima of the potential for 




FIG. 5: Plots of the minima of the potential. The colors and 
shapes of the points are used to identify distinct values of 
the superpotential at the corresponding minima. Solid and 
dashed lines indicate BPS and non-BPS sectors, respectively. 
The dot-dashed lines represent BPS sectors whose orbits in 
field space go through infinity and are never realized in prac- 
tice. The drawings start at the left top, and are given for 
some tipical values of b, in the intervals 6>1,0<6<1, 
-1/3 < b < 0, -1 < b < -1/3, -3 < 6 < -1, 6 < -3 
respectively. 



several values of b. The colors and shapes are used to 
identify distinct values of the superpotential at the cor- 
responding minima. The horizontal and vertical axes rep- 
resent the (j) and x fields and we indicate the BPS and 
non-BPS sectors with solid and dashed lines, respectively. 
However, these lines should not be confused with orbits 
in field space which in general are curved. 

In Fig. 4 we represent the specific cases with 6 = 1, 
6 = 0, b = —1/3, b = —1 and b = —3, which are special 
cases either because one of the BPS tensions vanishes 
or correspond to a case where the Z4 symmetry is re- 
stored {b = 1). We will see later that these values of b 
correspond to critical points where the system changes 
behavior. For 6 = 0, the horizontal red and green lines 
show the continuum of minima with distinct values of 
the superpotential. Note that, for 6 = 1, the potential is 



FIG. 6: Plots of the BPS orbits for b = 0.5 (left panel) and 
b — —0.3 (right panel). Notice that almost all the orbits are 
curved. The figure for 6 = —0.3 is illustrative of the fact that 
for b < there are no finite paths in field space connecting 
the minima of the BPS sectors represented by a dot-dashed 
line in Fig. 5. 






p 




FIG. 7: Plots of the final results of numerical simulations for 
6 = 1.2 (top right), b = 0.8 (middle left), b = -1/3 (middle 
right), 6 = — 1 (bottom left) and 6 = — 3 (bottom right), 
starting from the configuration shown in the top left panel. 
The initial configuration for all these simulations is shown 
in the top left panel. The minima, vi, V2, vs and V4, are 
represented respectively by the colors dark blue, light blue, 
green and orange (black, dark grey, light grey and grey in 
black and white). 



invariant under the transformation (0, x) (x? 0)- How- 
ever, the superpotential (22) is not invariant under this 



transformation (see top left panel of Fig. 4). 
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FIG. 8: Domain wall network simulations for 6 = 0.98. At 
the initial time the minima vi and vs were randomly assigned 
with equal probability to the grid points. The time deriva- 
tives of the fields, and were initially set to zero. We see 
that, despite not being present in the initial configuration, 
the minimum V4 dominates at late times leading to a natural 
solution to the cosmological domain wall problem. 



In Fig. 5 we represent some generic cases, for b at some 
tipical values inside the several regions delimited by the 
values shown in Fig. 4. The dot-dashed lines represent 
BPS sectors whose orbits in field space go through infinity 
and are never realized in practice (these orbits must pass 
through a point where the energy density diverges) . The 
BPS orbits for b = 0.5 and b = —0.3 are given in Fig. 6. 
(left and right panels respectively). Notice that almost 
every BPS orbit is curved. For b = 0.5, and in fact for 
any 6 > 0, there are BPS orbits connecting every minima, 
except Vi and ^3. On the other hand the right panel of 
Fig. 6 is illustrative of the fact that, for 6 < 0, there 
are no finite paths in field space connecting the minima 
of the BPS sectors represented by a dot-dashed line in 
Fig. 5. 

We now study bifurcation and the possibility of pat- 
tern changing as we vary the parameter b of our model. 
To do this, we numerically simulate the model for some 
values of 6, namely b = 1.2, b = 0.8, b = —1/3, b = — 1 
and b = —3. All the simulations performed in this paper 
are 1024^ two-dimensional simulations in a matter domi- 
nated universe. The expansion of the universe introduces 
a damping term which slows down the domain walls and 
damps the radiation produced by them. The initial con- 
figuration for all these simulations is shown in Fig. 7 (top 
left panel). The simulation results for b = 1.2, b = 0.8, 
6 = — 1/3, 6 = — 1 and b = —3 are shown in Fig. 7 from 
left to right and top to bottom respectively, starting at 
the top right panel. Note that the results for b = 1.2, 
b = 0.8, b = —1/3 and b = —3 are representative of the 
behavior in the intervals b > 1,0<6< 1, —1 < b < 
and 6 < — 1 respectively. The small curvature of the do- 




FIG. 9: Domain wall network simulations for b = 0.98. Ini- 
tially, at each grid point, one of the minima {vi to V4) was 
randomly assigned with equal probability. The time deriva- 
tives of the fields, (f) and were initially set to zero . In this 
case the junctions are of the X-type but the minimum V4 still 
dominates at late times. 



main walls in the final configurations are indicative that, 
by the end of the simulations, the domain walls have 
slowed down considerably but are not yet frozen. Some 
very slow dynamics is still going on. 

For 6 > 0, X-type junctions are preferred. Notice that 
r24 = T12 + ri4 and consequently, in one dimension, the 
diagonal solution in the sector V2 ^ V4 has exactly the 
same energy as the sum of the two edge solutions in the 
sectors vi ^ V2 and vi ^ V4. Still, the discussion in 
Section II implies that X-type junctions are stable for any 
non-zero intersection angle. We have also numerically 
verified that, for any value of 6 7^ the direct path in the 
non-BPS diagonal sector vi ^ vs is always disfavored 
on energetic grounds. Hence, for 6 > only edge walls 
appear in the simulations. The intersection angles are 
related by ri2Cos(^4/2) = ri4Cos(^2/2) where O2 and O4 
are the angle subtended at the junction by the domains 
2 and 4, respectively (note that ri2 = r23 and ri4 = r^^). 
Consequently, the variation of the intersection angle, ^2, 
of the X-type junction with 6 > is a consequence of the 
dependence of the ratio between ri4 + rs4 and ri2 + T23 
on the parameter 6 (if 6 = 1 then O2 = O4). Notice (see 
Fig. 3, bottom panel) that the diagonal domain wall with 
tension ris decays into two domain walls with tensions 
ri4 and T34 {b < 1) or ri2 and r23 {b > 1). 

For b < there are five BPS-sectors but for two of 
them (those represented by dot-dashed lines in Fig. 5) the 
corresponding domain walls are never realized in practice. 
Consequently, if 6 < only the domain walls connecting 
the sectors represented by a solid line in Fig. 5 appear in 
the simulations, as it is clearly shown in Fig. 7. The case 
with b = —1 (bottom left) is special in the sense that the 
minima V2 and V4 coincide. 
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FIG. 10: Evolution of the walls comoving area as a function 
of the conformal time, 77, in 1024^ matter era simulations cor- 
responding to the two situations illustrated in Figs. 8 (dashed 
lines) and 9 (solid lines). In this case, the asymmetry param- 
eter is b = 0.70, 0.86, 0.92, 0.96 from left to right, respectively. 

C. The cosmological domain wall problem 

In the simplest scenarios domain walls are known to 
evolve towards a scaling regime, in which the character- 
istic scale of the network is of the order of the Hubble 
radius, until they dominate the energy density of the uni- 
verse. In fact, in order to avoid catastrophic cosmologi- 
cal consequences, in standard domain wall scenarios, the 
symmetry breaking phase transition that produces the 
domain walls has to be smaller than about 1 MeV. Bi- 
ased domain wall networks may provide a way to evade 
the Zel'dovich bound [18' by suppressing domain walls 
before they have any chance of playing a cosmological 
role. This may be achieved if the domain walls separate 
regions with different values of the vaccum energy den- 
sity [ini [20] , in which case the evolution tends to favor 
the lower energy minimum if the biasing is large enough. 
This is the mechanism behind the devaluation scenario 
[2T] which has been proposed as a possible solution to 
the cosmological constant problem and has recently been 
studied in detail in ref. [22]. Another possible solution to 
the cosmological domain wall problem is based on initial 
conditions which favor one minima over the others [19 . 

Here, we explicitly show that small asymmetries of the 
vacuum potential may lead to a suppression of domain 
walls connecting specific minima of the scalar field po- 
tential thus leading to a biasing mechanism responsible 
for the suppression of the domain wall network. In Fig. 
8 we plot the results of a network simulation for b = 0.98. 
At the initial time the minima vi and ^3 were randomly 
assigned with equal probability to the grid points and 
the time derivatives of the fields, (f) and were initially 
set to zero. The domain wall connecting the minima Vi 
and vs is disfavored and an asymmetry is created which 



favors the dominance of domain 'U4 over the others. In 
Fig. 9 we plot the results of a domain wall network sim- 
ulation of the same model but in this case all the four 
minima were initially assigned to the grid points with 
the same probability. The junctions are of the X-type 
but the minimum V4 still dominates at late times due to 
the asymmetry of the potential which leads to the differ- 
ent domain wall tensions. We have also verified that the 
results for b = 1.02 are very similar to those obtained for 
b = 0.98 except that, in this case, it is the minima V2 
(rather than V4) which dominates over the other three. 
Note that, in both cases we have a multi-tension network 
whose evolution tends to favor domain walls with lower 
tensions thus introducing a bias which favors domains 
surrounded by those walls. We numerically verified that, 
for 6 < 0, the domain wall network has no junctions and 
the evolution is rather trivial, since it tends to favor the 
only minimum of the potential which connects all other 
minima with BPS domain walls. Of course, this was al- 
ready expected taking into account the above discussion 
and the results shown in Fig. 7 for 6 < 0. 

In Fig. 10 we plot the evolution of the comoving area 
of the domain walls as a function of the conformal time, 
7^, corresponding to the 1024^ matter era simulations pre- 
sented in Figs. 8 and 9 (dashed and solid lines respec- 
tively). However, in this case, we have considered sev- 
eral different possibilities for the asymmetry parameter 
{b = 0.70,0.86,0.92,0.96 from left to right respectively). 
Fig. 10 clearly shows that the domain wall network is 
exponentially suppressed after some time. These results 
are similar to those obtained in Ref. [19 where other bi- 
asing mechanism were considered. The total suppression 
of the domain walls occurs for both sets of initial con- 
ditions. However, in the case where all the minima are 
initially randomly assigned to the grid (solid lines), this 
suppression happens for slightly larger values of the con- 
formal time than in the case where only the minima Vi 
and Vs are initially present (dashed lines). As expected, 
for larger values of |6 — 1| this suppression occurs at ear- 
lier times. In particular, b = 0.96 roughly corresponds 
to the smallest asymmetry of the potential for which the 
exponential suppression of the domain wall network still 
occurs on the trustfull part of the dynamical range of 
these simulations (where the particle horizon is smaller 
than half of the box size). See Refs. [12 for technical 
details of the simulations. 



IV. CONCLUSIONS 

In this work we have investigated specific models de- 
scribed by two real scalar fields and studied the role of 
bifurcation in changing domain wall network geometry 
and topology. First we have considered a class of mod- 
els without a superpotential controlled by a small pa- 
rameter, e. In this case a perturbative expansion up to 
first-order in e or the computation of the tensions asso- 
ciated with straight line orbits allows for the analytical 
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determination of the domain wall tensions in the |£| <C 1 
limit. We tested the analytical results by comparing with 
numerical simulations and found a very good agreement 
for \s\ < 0.2. We also confirmed that the model bifur- 
cates which, in more than one spatial dimension, results 
in a change from a phase where F-type junctions are fa- 
vored {e > 0) to another one with only X-type junctions 
{e < 0). 

Note that there are significant differences on the evo- 
lution of the macroscopic properties of domain wall net- 
works which depend on whether or not junctions are 
present and, in the former case, on the dimensionality of 
the junctions. Although the presence of junctions does 
not seem to help much from the point of view of making 
domain walls a viable dark energy candidate [12^, their 
impact on the dynamics of the network needs to be taken 
into account in any quantitative study of cosmological 
consequences of domain wall evolution. A good example 
are models in which the domain walls separate domains 
with different values of fundamental couplings. 

In this paper we have also considered generic models 
with two real scalar fields with a quartic potential gov- 
erned by a cubic superpotential. In this case we have 
shown that the number of parameters affecting the ra- 
tios between the various domain wall tensions can be 



reduced to 4. We then studied, in detail, a particular 
subset of models controlled by a single parameter which, 
despite its simplicity, has revealed a very rich structure, 
with several regions of distinct behavior. We have also 
shown that these models incorporate a natural solution 
to the cosmological domain wall problem. If the Z4 sym- 
metry of the scalar field potential is broken into Z2, by 
moving b away from unity, a bias is created which tends 
to favor one minimum of the potential over the others. 
This type of mechanism might operate in realistic domain 
wall models, leading to a suppression of the domain walls 
before they become cosmologically relevant. 
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